DESY 97-210 



ISSN 04118-9833 



Thrust Jet Analysis of 
Deep— Inelastic Large— Rapidity— Gap Events 



HI Collaboration 



Abstract 

A thrust analysis of Large-Rapidity-Gap events in deep-inelastic ep collisions is pre- 
sented, using data taken with the HI detector at HERA in 1994. The average thrust 
of the final states X, which emerge from the dissociation of virtual photons in the 
range 10 < < 100 GeV'^, grows with hadronic mass Mx and implies a dom- 
inant 2-jet topology. Thrust is found to decrease with growing Pt, the thrust jet 
momentum transverse to the photon-proton collision axis. Distributions of are 
consistent with being independent of Mx- They show a strong alignment of the 
thrust axis with the photon-proton collision axis, and have a large high— Pj tail. The 
correlation of thrust with Mx is similar to that in e~^e~ annihilation at \/See ~ Mx, 
but with lower values of thrust in the ep data. The data cannot be described by 
interpreting the dissociated system X as a qq state but inclusion of a substantial 
fraction of qqg parton configurations leads naturally to the observed properties. The 
soft colour exchange interaction model does not describe the data. 
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1 Introduction 



The Large-Rapidity-Gap (LRG) events observed at the ep coUider HERA ||T|, ^ are attributed 
mainly to diffraction Q]. The kinematics of these events can be discussed quite generally in 
terms of the process ep e'XY (see Fig. 1), where the two hadronic systems X and Y are 




Figure 1: Illustration of the process, ep e'XY, in which the hadron systems X and Y are separated 
by the largest rapidity gap in the event. W is the invariant mass of the colliding virtual photon - beam 
proton system. —Q^ is the square of the 4-momentum transfer at the (e, e') vertex, and t at the {p, Y) 
vertex. Mx and My are the masses of systems X and Y, respectively. 

separated by the largest rapidity gap in the event. In the HI experiment, events are selected 
with no activity over a particular large pseudorapidity range adjacent to the outgoing proton 
beam (forward) direction. In these events the system Y remains as a proton, or is excited to a 
low mass baryonic system by a colourless exchange with a small transverse component in the 
momentum transfer Q iP = p—Y. The incoming photon, ranging from very low to high virtuality 
Q^, can be excited ^ into a vector meson (7W +p^V + Y), or can dissociate into a high 
mass state X (7^*) +p^ X + Y). 

In this paper, the topological structure of final states, which emerge from the dissociation 
of highly virtual photons, is investigated. In the photon dissociation (PD) picture this process 
can be related to the fluctuation of the virtual photon into partonic states |^ 0, followed 
by parton-proton quasi-elastic diffractive scattering, as shown in Fig. 2. In the proton infinite 
momentum frame (IMF) picture the diagrams in Fig. 2 can also be described in terms of 
deep-inelastic scattering (DIS) off an exchanged object with partonic content |3|, |^. 

The simplest PD process was first discussed by Bjorken and Kogut @, in terms of the quark- 
parton model in the form of the Aligned Jet Model (AJM). Here the photon fiuctuates into a 
quark anti-quark pair long before it actually collides with the proton. 

Quark anti-quark (qq ) configurations with large parton momenta kf transverse to the colli- 
sion axis are predicted to have a reduced cross section in comparison with the expectation for 
a hypothetical collision of two single quarks with the proton. This suppression of large kt has 



been explained in terms of "colour transparency" |11], since the physical transverse separation 
of the qq state is small on average, leading to a small colour dipole moment, or equivalently 
to mutual screening. However, qq configurations with a small transverse momentum difference 



^Here, and throughout this paper, the term pomeron and the symbol IP are used for the colourless exchanged 
object independently of its nature, unless the difference between a pomeron and a subleading reggeon exchange 
is significant in the context. 

^The names for the respective processes as given in the caption are appropriate for the IMF interpretation of 
these diagrams. Because the two interpretations of the diagrams in Fig. 2 are equivalent, these names are used 
also when referring to each of Fig 2a through 2d in the PD interpretation. 
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Figure 2: Diagrams of parton processes for LRG events: Born term diagram (a), boson-gluon fusion 
(BGF) (b), QCD-Compton (QCD-C) (c, d). The q ^ q exchanged configurations are impUed. 



do not lead to such screening effects and hence are not suppressed in their interaction with the 
proton. In the rest frame of X, the AJM picture thus leads to a back-to-back 2-jet configuration 
with small jet transverse momenta Pt relative to the incoming proton beam direction, typical of 
soft interactions. 

The QCD extension of this simple AJM model with gluons radiated by the qq system was 



proposed in |12], and has been studied further by several authors [0], in particular by invoking a 
soft colour interaction (SCI) process In the IMF view this corresponds to a QCD evolution 

of the structure function of the diffractive exchange. 

Two further and completely different approaches have been used to explain LRG event pro- 
duction. In one of them the LRG events are the result of normal DIS scattering with a subsequent 



SCI process |13], generating a large rapidity gap. In the other one it is proposed that the virtual 
photon fluctuates according to the Generalized Vector Meson Dominance Model (GVDM) [l 
into virtual vector mesons which then undergo diffractive scattering off the proton ||l5|, ^ 

All the above approaches are phenomenological and one can alter parameters to describe the 
inclusive production cross section for LRG events. They lead, however, to distinctive features 
in the topology of the hadronic final state. 

Topological studies were instrumental in analysing the transition of timelike photons into 
partons in e^e~ annihilation, even at energies y/See which were too low to resolve 3-jet from 
2-jet events [^]. The contribution of higher order parton configurations like qqg , relative to 
qq , was determined without applying jet algorithms, by measuring the deviation from colinear 
momentum flow using the event-shape variable thrust T |l8| , p!9[ . Data on the evolution of thrust 
with ^See also allowed a separation of hadronization from perturbative QCD effects to be made. 

In this analysis, thrust is used to study topological features of the hadronic final state in 
7* dissociation in a kinematic regime where diffractive contributions are dominant [^] . Here the 
final state partonic structure and the related hadronization effects may be different from those 
in e~^e~ annihilation events. The correlation of thrust with Mx may be exploited to separate 



these two contributions. In addition, the ahgnment, i.e. the correlation of the thrust direction 
with the incoming proton direction, can be investigated by studying the Pt distribution of thrust 
jets with respect to the proton direction and can be used as a further tool for the analysis of 
the dynamics of LRG events. 

Because of a forward rapidity gap requirement, the hadronic final states of accepted events 
are well contained in the central detectors of HI. This provides measurement conditions similar 
to those in e~^e~ annihilation experiments at symmetric e~^e~ colliders. No jet energy threshold 
is required in the thrust analysis and all events in the Mx range under study are included. An 



event shape analysis of LRG events has been reported recently by the ZEUS collaboration |21] 



2 Detector description 

In the following, only detector components relevant for this analysis are reviewed. A detailed 
description of the HI detector can be found elsewhere [p^ ]. 

The "backward" electromagnetic calorimeter (BEMC) has full azimuthal coverage, and ex- 
tends over the range 151° < 9 < 176°, where 6 is the polar angle with respect to the proton beam 
direction, as seen from the nominal beam collision point. The BEMC was used to trigger on 
and measure the energy of the scattered electron in DIS processes. The electromagnetic energy 
resolution is aE/E = 0.10/^^ [GeVJ © 0A2/E[ GeV ] © 0.03 while the absolute electro- 

magnetic energy scale is known with an accuracy of 1%. The BEMC hadronic energy scale is 
known to a precision of 20%. The backward proportional chamber (BPC), located in front of 
the BEMC detector, has an acceptance of 155.5° < 9 < 174.5° and, in conjunction with the 
interaction vertex, measures the direction of the scattered electron with a precision of 1 mrad. 
These "backward" detectors accept DIS processes with values ranging from 5 to 120 GeV^. 
The liquid argon (LAr) calorimeter extends over the polar angular range 4° < ^ < 154° with full 
azimuthal coverage. The electromagnetic energy resolution is cte/E ~ 0.11/ ^/E\GeV] © 0.01, 
while the hadronic energy resolution is cte/E « 0.50/ ■\/ E [GeV] © 0.02 as determined in test 
beams. A study of the transverse momentum balance between the hadronic final state and the 
scattered electron, Pt^h ~ Pte ■• shown that the absolute hadronic energy scale is known to 
an accuracy of 4%. 

The LAr calorimeter is surrounded by a superconducting solenoid providing a uniform mag- 
netic field of 1.15 T parallel to the beam axis in the tracking region. The track reconstruction 
is based on the information from the central jet chamber (CJC), the z-drift chambers and the 
forward tracker. These detectors cover a polar angular range of 5° < < 155°. 

Forward energy deposits at small angles are observed in several detectors near the outgoing 
proton beam direction. Particles reach these detectors both directly from the interaction point, 
and indirectly as a result of secondary scattering with the beam pipe wall or adjacent material 
such as collimators. These detectors are thus sensitive to particles well outside their nominal 
geometrical acceptances. The liquid argon calorimeter is sensitive to particles with pseudora- 
pidities r] = — In tan 0/2 up to r/ ~ 5.5. A copper-silicon sandwich calorimeter (PLUG) allows 
energy measurements to be made over the range 3.5 < < 5.5. The three double layers of drift 
chambers of the forward muon detector (FMD) are sensitive to particles produced at pseudora- 
pidities 5.0 < ry < 6.5. The proton remnant tagger (PRT), consisting of seven double layers of 
lead/scintillators, located 24 m from the interaction point, covers the region 6.0 < i] < 7.5. 

These detectors overlap considerably in their rapidity coverage, thereby allowing intercali- 
bration of their efficiencies. 



3 Physics models for LRG event production 



In the PD picture the photon fluctuates into a gg, qqg, or states with more partons, where at 
least one of the partons is virtual. A virtual parton, if it has a low kt relative to the 7*p collision 
axis, can scatter difFractively (quasi-elastically) with the proton yielding a real parton. In the 
qq configuration (Fig. 2a) this gives an aligned 2-jet event, and at least one aligned jet in the 3- 
parton configurations (Figs. 2b, c and d). It is expected qualitatively that the qqg configuration 
has a larger cross section than the qq configuration for difFractive scattering off the proton 

In the IMF picture the diagrams in Fig. 2 are interpreted as the DIS probing a colourless 
exchange object. This is implemented in the RAPGAP ||2^ Monte Carlo (MC) program. The 
RAPGAP model employs deep-inelastic electron scattering off pomerons [|lO| and off (subleading) 
reggeons coupling to the proton. The pomeron is ascribed a quark and a gluon content. The 
ratios of the BGF to the Born term and to the QCD-C contributions (see Fig. 2) depend on the 
gluon and quark contents of the exchanged objects. These have been determined from a QCD 
analysis with DGLAP evolution p^] of the diffractive structure function measured by the HI 



collaboration |2£], with the result that most of the pomeron momentum is carried by gluons. 



For the reggeon, the quark and gluon content is taken to be that of the pion p^]. 

The RAPGAP model is used in two modes. The matrix element mode (ME) includes the 
BGF and the QCD-C processes as first order perturbative QCD (pQCD) matrix elements, with 
propagators for the parton mediating the hard process, while the diagram in Fig. 2a (Dia- 
gram 2a), is treated as DIS off a quark in the pomeron. The ARIADNE |2J] mode (AR) also 



treats the BGF process as a matrix element, but treats Diagrams 2c and 2d as colour dipole 
radiation from the qq pair of Diagram 2a. A cut-off in terms of the square of the transverse 
momentum of the hard sub-process, > Ptmini— ^ GeV^) is applied to keep the contribu- 
tion from diagrams treated as pQCD matrix elements below the measured total cross section at 



each {Mx,Q'^,x) point. With this choice |24] of cut-off, the corresponding diagrams of Fig. 2 
contribute about 50% of the total diffractive cross section within RAPGAP. The remainder is 
generated with the process of Diagram 2a. The cut-off value can be raised, but it cannot be 
lowered without the need to readjust other parameters influencing the total cross section. The 
topological properties, like thrust or thrust jet Pt, of the simulated events are sensitive to the 
fraction of matrix element events and depend on the value chosen for pt^min- 

In the RAPGAP ME mode higher order QCD effects are included by means of leading- 
log parton showers for all final state partons. Predictions of this mode are labelled as 'RG 
ME-I-PS' in the following. In the AR mode higher orders are included as additional gluon 
radiation according to the Colour Dipole Model [^] (labelled as 'RG AR+CD'). In both modes 
the hadronization is as given by the Lund string fragmentation model 



In generating events corresponding to Diagram 2a the qq system is taken to be aligned 
exactly with the 7*iP axis, i.e. without intrinsic kf. The virtuality of the quark mediating the 
interaction between the photon and the pomeron is thus always taken as the smallest value 
which is kinematically allowed. These events can be considered as representing the simple AJM, 
which differs from the RAPGAP simulation of Diagram 2a by having no QCD dipole radiation, 
and by neglecting transverse momentum transfer to the proton. The corresponding predictions 
are labelled as 'RG qq'. 

The non-diffractive ep DIS plus SCI picture is implemented in the LEPTO event genera- 



tor |31]. This employs the same hard processes as shown in Fig. 2, but now off the proton. The 
ratio of the generated qqg to qq events is determined by an infrared cut-off ^ and the parton 
density functions for the proton. Both are adjusted to describe the total DIS event sample, but 



The infrared safety scheme in LEPTO avoids a sharp cut-off in pt- 



there is no parameter in LEPTO to adjust specifically for the topological properties of the LRG 
events. This is different from the RAPGAP model where the relative gluon and quark content 
of the pomeron and the pt,min value can be adjusted separately to describe the LRG event data. 
The parameter which controls the strength of the SCI process has no influence on the topological 
properties of the system X. The SCI process, and hence the LRG which separates the systems 
X and Y in Fig. 1, is not limited to occur between the remnant partons of the struck proton and 
those emerging from the hard process. Thus, the final state X can in general not be identified 
with a specific hard process. 

Predictions from the GVDM picture for final state properties require further specifications 
of the diffraction process as a hadronic collision. 



4 Thrust jet analysis method 

In the centre of mass of a system X of particles, the thrust method determines the direction of 
the unit vector a along which the projected momentum flow is maximal 111]. Thrust is computed 



as 119]: 



N N 



T = (1/^ \pi\) ■ maxJ2\Pi ■ «l 

i=i °- i=i 

where pi represents the momentum of particle i, in the rest frame of the A^ particles 0. 

Given the thrust axis a, the directional sense of which is arbitrary, the N particles can then 
be grouped uniquely into two subsets (thrust jets), depending on whether they belong to the 

hemisphere with positive or negative momentum component along the thrust axis. The summed 

Nj 

particle momenta of hemisphere / form the jet 4- momentum Pi = Pj with Nj the number of 

i=i 

particles in hemisphere /, for / = 1,2. The two thrust jets have independent masses, and equal 
but opposite 3-momenta: |-Pi,2| = P- 

Thrust values range from a maximum value of T = 1 in the case of a 2-particle state or any 
colinear configuration, to a minimum value of T = 0.5 obtained in an isotropic system X with 
infinite multiplicity. A symmetric 3-particle configuration yields a value T = 2/3 and leaves the 
direction a arbitrary in the 3-particle plane, while a non-symmetric topology gives T > 2/3 and 
a thrust axis pointing in the direction of the most energetic particle. A non-symmetric 3-particle 
topology will thus appear as a 2-jet like configuration with T < 1. 

In a multihadron state emerging from a partonic process, the two back-to-back reconstructed 
thrust jets are correlated with the hard partons in the following way. In the case of a 2-parton 
configuration, the thrust value as determined from the final state hadrons - i.e. at the hadron 
level - is smaller than 1, and the direction of the thrust axis remains parallel to the direction of 
the two partons to the extent that the hadrons can be correctly assigned to 'parent-partons'. For 
an underlying 3-parton system, the reconstructed thrust axis at the hadron level is correlated 
with the direction of the most energetic parton. This property has been verified to persist down 



* There is another definition of thrust (Tm) in the literature jl^, where the normalization is not to ^ \pi\ as 

j=i 

in T but rather to Mx, where Mx is the invariant mass of the system X. If final states contain massive particles, 
Tm < T. To compute Tm accurately, it is necessary to identify all final state particles, which is usually not 
possible experimentally. 



to final state masses Mx ^ 5 GeV, using qqg events from the RAPGAP event generator |24] 
with hadronization and detector effects included in the event simulation. 



A transverse thrust jet momentum Pt can be defined relative to a reference direction r as 
Pt = P ■ sin Q with the angle between Pi and r. In this analysis, the proton beam direction 
transformed into the centre of mass frame of the system X is chosen for r . Since the 4-momentum 
transfer squared t is small in LRG events, this direction is a good approximation to the 7*iP 
axis. For t = tmin they are identical. Monte Carlo studies have shown that a diffractive t — tmin 
distribution with an exponential slope parameter 6 = 6 GeV~'^ leads to an average smearing of 
Pt of less than 0.3 GeV. Unlike the 7* direction, when calculated from the scattered electron, 
the proton direction in the X rest frame is not affected by QED-radiation. The JP direction is 
not well determined experimentally. 

If events with large thrust jet Pt are dominantly caused by single gluon radiation from an 
aligned qq state, then large Pt will be correlated with lower average thrust (T). If the larger 
Pt values were instead caused by a rotation of a qq configuration to larger angles, i.e. by larger 
intrinsic kt, then thrust would not change with Pt. Hence, the observation of a decrease of (T) 
with increasing Pt is evidence for the presence of qqg contributions in the final state. 



5 Event selection and final state reconstruction 

The DIS event sample used in this analysis was obtained from data collected during the 1994 
HERA operation, and corresponds to an integrated luminosity of about 2pb~^. The event 
trigger required a minimum energy deposition in the BEMC of 4 GeV. The events were selected 
according to the kinematical properties of the scattered electron f\ [^2|| and the hadronic final 
state. The following variables are used: 

Q2 = 4EeE^, cos\ej2) ye = l- E^,/Ee sin2(0^,/2) x = QVives), 

calculated from the electron beam energy Ee{= 27.5 GeV), the energy E^' and the polar angle 
9^' of the scattered electron, and y/s, the invariant ep mass. The invariant mass Mx of the 
hadronic system X is calculated as a geometric mean : 

Mj, = Mx{e) Mx{h) (1) 

with 

Ml(e) = [E, - P,,, - {E^, - P^^MEh + P.,h) - PIh " 



Ml{h) = El-Pl^-Pl^-Pl^ 

where 

Ph 'Y^hadrons P'x{y,z),h ^hadrons P'x{y,z),i 

^ The data analysed here were taken with a positron beam in HERA. The term electron is generic. 



and Ei, Px,i, Py,i,Pz,i are the individual energy and momentum components of the hadronic final 
state in the laboratory system. 

The hadronic 4-momenta were calculated using the interaction vertex and the energy deposit 
clusters in the calorimeters. In addition the measured momentum from each track was included 
up to a maximum of 350 MeV/c to compensate for losses in the calorimetric measurement of 
charged hadrons. This limit was determined by the requirement of a balance between -P^g', 
the transverse momentum of the scattered electron, and Pt^h, the transverse momentum of the 
observed hadrons. For the LRG events, which have only very small unobserved pt contributions 
in the forward direction, this procedure yields a width of 1.2 GeV in the Pt^h — Pt e' distribution. 

The reconstruction errors for Mx(e) and Mx{h) are mostly independent. Large errors can 
occur because of initial state QED radiation [Mx(e)] or from acceptance losses in the backward 
direction [Mx(/i)]. Studies with simulated events show that using the geometric mean [equ.(l)] 
reduces extreme deviations from the true Mx- 

From the above quantities and with q the 4-momentum of 7* the variables (3 and xjp can be 
calculated as follows 

^ 2q-{p-Y) Q'^+Mj.-t 
_ q-{p-Y) _ Q'^ + Ml-t _ x_ 

In the "proton infinite momentum frame" , the variable x jp can be interpreted as the fraction of 
the proton momentum transferred to the photon by the exchanged IP, while /? is the fraction of 
the IP momentum carried by the parton coupling to the photon. 

The 4-momentum transfer squared, t, between the incident proton and the final state Y is 
small in LRG events since the forward detector selection amounts to |t| < 1 GeV^. Therefore t, 
which is not measured, can be neglected in the above formulae. 

The following selection criteria were applied to the scattered electrons of the LRG sample: 



1. The electron energy was constrained to be large {E^i > 14.4 GeV) by requiring ye < 0.5. 
This keeps the photoproduction background at a negligible level, and removes events where 
the hadronic system X is boosted strongly towards the BEMC and the backward beam 
pipe hole. 

2. The electron scattering angle was limited to 156° < < 173° to avoid events with electron 
energy deposition near the edges of the BEMC. 

3. was required to be in the range 10 to 100 GeV^. In this range there is almost 
constant acceptance for electrons. 



Further selections were made with quantities derived from the full final state: 



1. The interaction vertex was required to be within ±30 cm of its nominal position and to 
be reconstructed with at least 1 track in the central jet chamber. 

2. xjp < 0.05 was required to suppress non-diffractive contributions. 

3. The largest rapidity gap was required to include the pseudorapidity range of the forward 
detectors (see section 2) by demanding the absence of significant activity in these detectors 
and, for rj > 3.2, in the LAr calorimeter. These requirements have a high efficiency to 
reject events with particles in the region 7.5 > rj > 3.2 |3S]. 



4. The mass Mx of the hadronic final state was required to be bigger than 4 GeV, the 
approximate hmit below which the thrust axis ceases to be correlated with a 2-parton 
axis. 

5. Mx < 36 GeV was required to avoid events where the rapidity gap requirement implies 
an acceptance below 25%. 

With these selection criteria 6865 events were collected in the Mx range from 4 to 36 GeV. 
They are distributed over the Mx intervals as shown in Table |l[ These events cover the range 
10-4 < X < 3.2 • 10-2, and 30 < W < 210 GeV. 





Mx intervals (GeV) 


4-6 


6-8 


8-11 


11-15 


15-19 


19-24 


24-36 


No. of events 


1439 


1278 


1404 


1130 


711 


502 


401 



Table 1: Number of observed events per measured Mx interval. 



6 Corrections for resolution and acceptance 

Using simulated events, the measured distributions have been corrected for detector resolution 
and acceptance losses to give cross sections in a kinematic range defined by 

My < 1.6 GeV^ xjp < 0.05, |t| < 1 GeV^, 10 < < 100 GeV^ and y < 0.5, (2) 

for the same Mx intervals as given in Table [l|. Events with these My , x jp and t limits have large 
rapidity gaps in the detector for kinematic reasons This cross section definition only allows 
for small non-diffractive contributions, which are typically less than 10% [^], and it avoids the 
need for a statistical subtraction from the data of a model-dependent non-diffractive background. 
It also reduces the systematic error due to model dependence of acceptance corrections. 

Two different Monte Carlo event simulations were used to correct the observed experimental 
distributions. These were RAPGAP in the ME-I-PS and in the AR-I-CD mode, but with simple 
(flat) quark and gluon distribution functions for the pomeron and without evolution. The 
generated events were passed through a full detector simulation and had the same event recon- 
struction as real data. The Monte Carlo simulation of this simplified model described quite well 
the kinematic distributions of Q^,xip and W, and further variables describing the detector re- 
sponse. The observed Pt distributions were described qualitatively. With an additional smooth 
weighting function for the generated events, which depends on thrust jet Pt at the hadron level, 
the Pt distributions in all Mx intervals were very well approximated over the full Pt range. The 
measured thrust distribution is also well reproduced after the weighting in Pt, as shown in Fig. 3. 

The contents of the 2-dimensional bins of the experimental {Pt, Mx) and (T, Mx) distribu- 
tions were multiplied by the ratios of the numbers of generated to reconstructed and selected 
events in the corresponding intervals for the MC distributions. 

Monte Carlo studies show that the thrust jet Pt is resolved to better than 1 GeV. The 
correction factors vary only slowly from bin to bin in all Pt distributions and are large (> 2) 
only for masses Mx > 15 GeV and for P^ < 1 GeV^. In the thrust distributions, the correction 
factors in each mass bin vary less than a factor two, except for the lowest and the highest thrust 
bins. Since those intervals contain only a very small fraction of the events (see Fig. 3), these 
factors have very little influence on the average thrust value. The corrected average thrust is 
calculated from the corrected thrust distribution. 

To estimate the systematic errors, the following contributions were taken into account: 
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Figure 3: Normalized, uncorrected thrust distribution, with the number of observed events, 
for data (fuU points) and the simphfied, Pt-weighted MC simulation (open points), as described in 
the text. 



uncertainty in the LAr hadronic energy scale; 

uncertainty in the BEMC hadronic and electromagnetic energy scales; 
uncertainty in the correction factors due to freedom in the Pt weighting; 
differences between event generation schemes (parton shower, colour dipole). 



The two latter contributions were the largest sources of systematic error. Their effect was 
determined for each data point from the differences between the four correction factors coming 
from the two event generation schemes - with colour dipole radiation or with parton showers - 
and each of them in the unweighted and the weighted mode. The unweighted Pt distributions 
at detector level of the simulated events were slightly flatter than the data, and the weighted 
ones were slightly steeper. 



Radiative corrections have been applied using the programs HERACLES [34| and RAPGAP 
in the matrix element and ARIADNE modes, as described in Section 3. The corrections are 
smaller than the statistical errors for all data points shown. 



The statistical, and the statistical and systematic errors added in quadrature, are indicated 
in the figures as the inner and the outer error bars, respectively. Systematic and statistical 
errors are given separately in the respective tables. 



7 Results 



7.1 Correlation of thrust with Mx,Pt and 



In Fig. 4a thrust is shown as a function of 1/Mx, averaged in seven Mx intervals (see Table 
in the range 4 to 36 GeV, and averaged over the kinematic region of equ.(2). The same data are 
also given in Table 0. The correlation of thrust with 1/Mx rather than Mx is shown because it 
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Figure 4a: (a) Average thrust (T), as a function of 1/Mx for LRG data, and for e+e~ 
annihilation data ^ as a function of l/y/See] 4b: Average thrust (T), as a function of 1/Mx, 
for all events, for events with Pt < 1 GeV and for events with Pt > 1 GeV, with predictions from 
RAPGAP (see text). The data points are given at l/{Mx). 



(Mx)(GeV) 


(T) 


stat. error 


syst. error 


B 


4.97 


0.777 


0.004 


0.008 


0.18 


6.98 


0.781 


0.004 


0.005 


0.16 


9.40 


0.802 


0.004 


0.006 


0.18 


12.81 


0.822 


0.004 


0.005 


0.16 


16.82 


0.847 


0.005 


0.006 


0.12 


21.20 


0.865 


0.005 


0.008 


0.10 


28.64 


0.875 


0.006 


0.009 


0.10 



Table 2: Mean thrust values as a function of (Mx)- The following bin limits for Mx were used: 4, 6, 8, 
11, 15, 19, 24, 36 GeV. To combine values of (T) from different Mx intervals, the relative cross section 
weights B are to be used. 



allows perturbative and hadronization effects on thrust to be distinguished in a simple manner. 
For Mx oo the latter should become negligible. For comparison, Fig. 4a includes the available 
data for thrust in e^e~ annihilation events 



17, 3q] plotted at Mx = -y/s^- The following features 



are noteworthy: 

1. (T) increases with Mx- It is also observed in these data that the particle multiplicity (not 
shown here) also increases with Mx- This implies that the final state is not isotropic, since in that 
case (T) would fall, and approach 0.5, as 1/Mx — >■ 0. The growth of (T) signifies an increasing 
back-to-back correlation of the momentum flow as the final state mass increases. Hence, the 
thrust direction corresponds to an effective topological axis of the event. This feature rules out 
fireball type models [^6| with isotropic configurations of final state hadrons in virtual photon 
diffractive dissociation. Similar effects have also been observed in hadronic single diffractive 



dissociation |l37 |. 

2. There is a striking similarity between the ^/s^ dependence of (T) in e^e~ annihilation 
and its Mx dependence in virtual photon dissociation, for the kinematical range of this analysis. 
The slopes are comparable, but (T) of DIS LRG events is lower by an approximately constant 



value of 0.025. This similarity is not so evident in |21]. 

For 4 < Mx < 36 GeV the thrust-mass correlation may be described by 

(T) = Too - H/Mx + F/Mi , (3) 

as first observed in e^e" annihilation experiments There Too is related to thrust at the 

parton level and hence to Og in the See range of the data. The coefficient H describes 
the contribution to thrust of jet broadening by hadronization. Such a power correction term 

confirmed by data , reproduced with simulated events using 
and is now discussed in the context of non-pertubative QCD |40| ]. 
Contrary to the situation in e~^e~ annihilation, a theoretical justification of the Too and the 
1/Mx terms has not yet been given for LRG events, but the data suggest that equ.(3) applies 
here also. The 1/M'^ term is an ansatz to account for a bias of thrust towards large values when 
such a maximum search method is applied to low particle multiplicity final states which occur 
at low Mx ■ In a fit of this parameterization to the Mx range of the LRG data the parameter 
Too is correlated strongly with H and F. To express the measured (T) dependence on Mx with 
minimally correlated parameters an expansion of (T) in 1/Mx at a mass Mx = Mq inside the 
data range is appropriate: 



was theoretically predicted 
limited pt hadronization 




(T) 

Here Tq is (T) at Mx = 



n-H- (1/Mx - 1/Mo) + F ■ {I /Ml - I /Ml) 

20 GeV, the coefficient F for the e' 



(4) 



Mq. In a fit, with M, 







Data 


To at Mx = 20 GeV 


H (GeV) 


ep DIS LRG 
e^e~ annihilation 


0.857 ± 0.003 
0.881 ± 0.002 


1.78 ±0.03 
1.72 ±0.05 



Table 3: Results from the T — Mx correlation (eqn. 2) fitted to this data and to e+e data |3^, in 
the range from 6 < Mx < 36 GeV. Tq is the value of (T) at Mx = 20 GeV, and H the coefficient of the 
hadronization (l/Mx) term. The errors given are statistical errors for the ep data and statistical and 
systematic errors combined for the e'^e" data. For the ep data Tq has in addition a systematic error of 
±0.006. The H values for ep and e~^e~ have in addition a common error of 0.11 from their correlation 
to the uncertainty in the low multiplicity fluctuation term F. 

annihilation data was poorly constrained and was thus fixed to the result found in LRG events 
(F = 5.0 ± 0.6 GeV^). The results of the fits in Table | show that the power corrections {H) 
are compatible, while the thrust values are significantly lower in LRG than in e~^e~ events. 

The result {T) ^j^q < (T)^^ means that the final state in DIS LRG events cannot be under- 
stood as a simple qq state with standard limited pt hadronization |3£]. It implies that higher 
parton multiplicities are even more important in DIS LRG events than in e~^e~ annihilation. 
For a qq parton configuration, {T)j^j^q can only be less than (T)gg if the power corrections are 



stronger than in e"'"e~ annihilation data. Then the value of H would have to be larger than found 
in this fit. A quantitative interpretation of the difference {T)j^j^q — {T)^^ in the PD picture in 
terms of contributions from the Born-term, and from the QCD-Compton and BGF subprocesses 
(see Fig. 2), requires a treatment of the non-perturbative diffractive vertex in these diagrams. 
This is not yet available. In the IMF picture, the expressions for {T) j^^q are divergent. 

Contributions from underlying configurations with three, as in Fig. 2b, c and d, or more 
partons can also be studied using the correlation between Pt and thrust. In Fig. 4b we show (T) 
for events with Pt smaller and larger than 1 GeV. At each Mx point average thrust for large 
Pt is lower than for small Pt , demonstrating the presence of configurations with at least 3 (non- 
colinear) partons. The observed difference of (T) cannot be attributed to the selection of events 
with different hadronization properties. If this were the case, the difference would have to be 
proportional to 1/Mx, and this is not compatible with the data. On the contrary, the difference 
is compatible with a higher, Mx independent, parton multiplicity. Monte Carlo studies confirm 
that this Pt dependent (T) behaviour is that expected with 3-parton configurations, as seen in 
the curves in Fig. 4b, which show RAPGAP (ME+PS) expectations for a mixture of qq and 
qqg configurations. 

The dependence of (T) on was also studied. The value of (T) changes by no more than 
±0.01 over the range 10 to 100 GeV^. Since {Mx) is found to be independent of Q^, it can 
be concluded that the data show almost no correlation between thrust and Q^. 

7.2 distributions 

The cross sections differential in Pt are shown in Fig. 5a and in Table ^ for six Mx intervals in 
the mass range from 6 to 36 GeV, and normalized to unity in each Mx interval. The Mx interval 
from 4 to 6 GeV is not included since the flattening of (T) with 1/Mx in Fig. 4a, indicates that 
effects from low multiplicity fluctuations are becoming strong. For each Mx interval the data 
are displayed only in those bins of Pt that are kinematically fully accessible. In Fig. 5b the Pt 
distribution of one mass range (19 < Mx < 24 GeV) is displayed on a linear scale in Pt - The 
data show: 

1. a steep rise of the distribution towards low Pt values, demonstrating a dominant alignment 
of the final state thrust axis with the initial 7*p collision direction; 

2. a shape of the Pt distributions nearly independent of Mx, suggesting the factorizable form 

daiMx,P?)/dP? = AiMx) ■ BiP,^); (5) 

3. a substantial high Pt tail. 

The large high-P^ tail supports the conclusion in Sect. 7.1 that there are significant contri- 
butions from parton multiplicities higher than qq , such as those shown in Figs. 2b-d. The same 
conclusion may be gained from Fig. 6 where the fractions of events with Pt values above 1.0 GeV^ 
and above 3.0 GeV^, respectively, are displayed as a function of Mx- The fraction at the lowest 
Mx point in Fig. 6b is relatively low because of the reduced phase space for P^ > 3 GeV^ here 
(see Table Q). For higher masses the fractions reach about 64% for Pt > 1 GeV^ and 35% 
for Pt > 3 GeV^. The model predictions indicate that such large values of Pt require hard 
processes, which involve more than two final state partons, to play an important role in DIS 
LRG events. 
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Figure 5a: Normalized thrust jet distributions for six Mx intervals. The bin limits are 0, 
1, 3, 6, 12, 25, 43, 79, and 151 GeV^; 5b: Normalized thrust jet P^ distribution for 19 < Mx < 24 
GeV with a linear P^ scale. 
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Figure 6a: Fraction of events with P^ > 1 GeV^ for six Mx intervals, together with four model 
predictions (see text); 6b: as in a, but with P^ > 3 GeV^. 



8 Comparison with models 



The GVDM model in combination with quasi-elastic diffractive vector meson proton scatter- 
ing lie] leads to final states X with spin J = 1. For such a situation the width of the Pf 



{P'f) GeV^ 


iMx) (GeV) 


il/a){da/dP?) GeV-^ 


stat. error 


syst. error 


0.37 


6.96 


0.527 


0.018 


0.038 


1.88 


6.96 


0.144 


0.009 


0.022 


4.21 


7.10 


0.042 


0.004 


0.007 


0.38 


9.39 


0.470 


0.017 


0.039 


1.84 


9.36 


0.129 


0.008 


0.022 


4.30 


9.31 


0.052 


0.004 


0.007 


8.38 


9.54 


0.016 


0.002 


0.003 


0.40 


12.72 


0.451 


0.019 


0.033 


1.83 


12.74 


0.117 


0.009 


0.013 


4.22 


12.77 


0.045 


0.004 


0.005 


8.54 


12.78 


0.019 


0.002 


0.003 


15.97 


13.10 


0.0051 


0.0007 


0.0013 


0.41 


16.90 


0.404 


0.024 


0.038 


1.79 


16.78 


0.133 


0.012 


0.016 


4.32 


16.70 


0.053 


0.006 


0.008 


8.47 


16.98 


0.016 


0.002 


0.004 


17.10 


16.89 


0.005 


0.001 


0.001 


30.75 


16.85 


0.0007 


0.0002 


0.0003 


0.37 


21.01 


0.444 


0.030 


0.044 


1.84 


21.10 


0.111 


0.013 


0.016 


4.22 


21.40 


0.045 


0.006 


0.009 


8.08 


20.91 


0.020 


0.003 


0.005 


16.19 


21.07 


0.0038 


0.0008 


0.001 


33.18 


20.79 


0.0008 


0.0003 


0.0004 


51.70 


22.15 


0.00039 


0.00014 


0.00022 


0.41 


27.79 


0.361 


0.032 


0.039 


1.83 


28.05 


0.141 


0.018 


0.028 


4.26 


27.94 


0.048 


0.008 


0.013 


8.42 


28.26 


0.017 


0.003 


0.004 


17.38 


28.65 


0.0049 


0.0010 


0.0011 


32.09 


29.35 


0.0014 


0.0004 


0.0006 


58.28 


28.18 


0.00054 


0.00018 


0.00025 


113.20 


30.45 


0.00011 


0.00006 


0.00011 



Table 4: The {1 / cj){da / dP'^) distribution. The foUowing bin hmits were used: 0, 1, 3, 6, 12, 25, 43, 79, 
151 GeV^ for P^, and 6, 8, 11, 15, 19, 24, 36 GeV for Mx- The table gives the mean values of P^ and 
Mx in each 2-dimensional bin, the differential cross section in this bin divided by the total cross section 
for the respective Mx bin, and the statistical and the systematic errors in the same normalized units. 
The weight factors from Table |^ arc to be applied to combine distributions from different mass bins. 

distribution is proportional to M^- Within present statistics the data do not support such a 
contribution. Thus the spin parity structure of X is neither pure nor dominantly vector. 

The SCI model as presently implemented in LEPTO reproduces the Mx dependence of (T) 
as shovi^n in Fig. 7, but fails to describe the fractions of events with Pf larger than 1 or 3 
GeV^ (Fig- 6). This model as proposed has no means with which to adjust this prediction and 
therefore it cannot describe the LRG events exclusively. 

The PD picture is qualitatively consistent with two independent features of the data. Firstly, 
the strong, Mx independent, alignment peak seen in Fig. 5 is expected because of colour trans- 
parency. Secondly, the cross section ratio for qqg to qq configurations is expected to be larger 
than in e'^e~ annihilation, mainly because of the large contribution for the BGF process, again 
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Figure 7: Comparison of model predictions (see text) with the Mx dependence of 
(T) . 



as a consequence of colour transparency. This is supported by the observation of {T)^j^q < {T)^^ 
in Fig. 4a, and by the large tail in Fig. 5. The qualitative expectations are confirmed by 
calculations in limited regions of phase space . 

Calculations of the IMF picture in the RAPGAP implementation are shown in Figs. 6 and 
7. The AR+CD mode, which does not employ matrix elements for the QCD-C process, sys- 
tematically predicts thrust too high and fractions of events with large Pt too low. The ME+PS 
predictions are closer than those for the AR+CD mode, but the agreement with the data is not 
completely satisfactory. There are, however, several means to adjust the RAPGAP predictions 
to the data. One is a modification of the pt cut-off scheme. The present scheme, with an Mx 
independent pt cut-off value, is responsible for the counter-intuitive behaviour of the prediction 
for the Mx dependence of the large Pt event fractions in the ME+PS mode. In view of these 
model implementation problems, one cannot conclude that this approach is unable to describe 
the data. 

The RG qq curves, which approximate the simple AJM, indicate the extent to which the 
virtual photon dissociation process to qq only is unable to describe the data. 



9 Summary and conclusions 



A thrust analysis of DIS large rapidity gap (LRG) events, that is of events attributable to the 
process -f*P ^ ^Y, with Mx > 4 GeV, My < 1.6 GeV, xjp < 0.05 and 10 < < iqo GeV^ 
reveals: 



a dominant two-jet topology of the final state X; 



• a clear alignment of the thrust axis with the incoming proton direction in the X rest frame; 

• a dependence of the average thrust value (T) on Pt, the thrust jet momentum transverse 
to the incoming proton direction: (T) decreases when Pt increases; 

• a Pf spectrum independent of Mx', 

• a sizeable fraction of events with large Pf (> 3 GeV^); 

and in comparison with e^e~ : 

• the influence of hadronization on (T) is compatible with e^e~ ; 

• (T) is smaller than in e^e~ at Mx = y/s'e 



These features together cannot be explained by a pure qq configuration. It has been shown 
in a model independent way that a substantial contribution of qqg and higher multiplicity 
parton configurations is required in the final state X. This contribution is greater than in 
e^e~ annihilation data. 

There are two equivalent pictures for diffractive LRG event production. In the first the virtual 
photon fluctuates into multi-parton states together with quasi-elastic parton-proton diffraction. 
In the second the electron scatters off partons which constitute the diffractive exchange. The 
expectations of both pictures are broadly in agreement with the results of this analysis. Two 
other mechanisms for LRG event production, GVDM in the particular combination with quasi- 
elastic vector meson proton diffraction, and non-diffractive ep DIS plus subsequent soft colour 
interaction, do not explain all topological properties of the observed events. 

The abundance of parton multiplicities higher than qq implied by topological event properties 
is consistent with the large gluon content in diffractive exchange found in QCD-fits to the 
inclusive production cross section of DIS LRG events |20]. 
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